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ABSTRACT 


The  shoreward  and  seaward  propagating  tidal  wave  components  were  detennined  using 
four  co-located  pressure  and  velocity  sensors  longitudinally  deployed  in  Elkhorn  Slough, 
Monterey  Bay,  CA,  to  describe  tidal  wave  evolution  and  distortion.  Elkhom  Slough  is 
short,  reflective  (-100%)  estuary  consisting  of  a  narrow,  gently  sloping  main  channel  and 
vast  marsh  and  mud  flats  near  the  landward  boundary.  The  amplitude  reflection  for  the 
astronomical  tidal  constituent  is  -90%,  whereas  the  distortion,  described  by  summing  all 
non-astronomical  tidal  amplitudes,  reflection  is  large  (>125%)  for  all  stations  stating  that 
the  seaward  tidal  wave  is  more  distorted  than  the  shoreward  tidal  wave.  It  was  found  that 
the  reflective  time,  defined  as  the  time  it  takes  the  tidal  wave  to  propagate  to  the 
landward  boundary  and  back,  varies  as  function  of  tidal  elevation.  Surprisingly,  the 
reflective  time  increases  with  increasing  tidal  elevation  contrary  to  the  hypothesis  that  the 
shallow-water  tidal  wave  phase  speed  increases  with  increasing  tidal  elevation,  resulting 
in  a  reduced  time.  The  landward  end  of  the  slough  has  vast  low-lying  marshes  and  mud 
flats  that  are  inundated  during  higher  tidal  elevations,  causing  the  effective  (average) 
wave  depth  to  decrease  with  increasing  tide  elevation,  which  explains  why  the  tidal  wave 
phase  speed  decreases  instead  of  increasing.  The  elevation-dependent  reflective  time 
caused  by  a  time-varying  wave  phase  speed  modifies  the  shape  of  the  seaward  tidal  wave 
relative  to  the  shape  of  the  shoreward  tidal  wave,  which  explains  the  evolution  and 
distortion  of  a  tidal  wave  in  a  reflective  slough. 
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I.  MOTIVATION 


As  the  focus  of  U.S.  Navy  operations  has  shifted  from  the  blue  water  to  the  brown 
water  to  meet  the  changing  threats,  it  is  important  that  the  war  tighter  has  a  vast 
understanding  of  the  physical  processes  associated  with  all  environments  to  obtain  a 
competitive  advantage  against  the  adversary.  To  gain  a  greater  understanding  of  the 
littoral  environment,  the  Office  of  Naval  Research  (ONR)  has  funded  several,  large, 
cooperative,  research  projects  to  study  inlets  in  the  U.S..  Specifically,  in  May- June,  2012 
an  experiment  at  the  New  River  Inlet,  known  as  the  River  and  Inlet  Dynamics  experiment 
(RIVET)  was  conducted.  RIVET  2  is  currently  scheduled  for  the  Columbia  River,  in 
May,  2013.  The  New  River  Inlet  is  approximately  200  m  wide  with  a  very  large  back 
bay.  While  river  is  in  the  name  of  the  inlet,  there  is  no  river  forcing  the  inlet  is  actually 
tidally  forced.  The  Columbia  River  is  a  different  system  in  that  it  is  approximately  2  km 
wide  and  has  strong  river  forcing.  A  third  system  which  is  different  from  the  systems 
described  above  is  Elkhorn  Slough  which  is  a  short  (10  km),  shallow  (3  m),  highly 
reflective  system,  where  the  tidal  wave  is  a  standing  wave.  From  a  short  experiment 
conducted  in  the  NPS  field  studies  (OC4210)  in  2009,  differences  were  noticed  in  wave 
characteristics  between  the  shoreward  and  seaward  tidal  wave.  Due  to  a  lack  of  statistical 
confidence  in  the  measurements  and  the  increased  desire  by  the  Navy  to  gain  a  better 
understanding  of  the  physical  environment  of  the  littorals,  a  small  field  experiment  was 
conducted  to  explore  the  tidal  wave  characteristics  of  a  highly  reflective  estuary. 
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II.  INTRODUCTION 


As  a  tidal  wave  propagates  from  the  open  ocean  through  the  inlet  and  upstream 
into  the  estuary,  it  distorts  due  to  changes  in  water  depth,  channel  geometry  and  frictional 
effects  (Friedrichs  and  Aubrey  1988;  van  Rijn  2010).  The  tidal  distortion  is  generally 
represented  by  overtides  (tidal  harmonics  of  astronomical  constituents)  and  compound 
tides  (nonlinear  sums  and  differences  of  all  tidal  constituents)  that  the  astronomical  tides 
develop  in  shallow  water  (Dronkers  1964;  Uncles  1981;  Speer  and  Aubrey  1985; 
Friedrichs  and  Aubrey  1988;  Nidzieko  2010).  The  tidal  wave  distortion  is  described  by 
the  tidal  asymmetry,  which  is  the  inequality  of  the  ebb  and  flood  tidal  currents.  The 
inequality  detennines  the  net  sediment  transport  and  flushing  patterns  (Dronkers,  1986; 
Friedrichs  and  Aubrey,  1988).  For  example,  an  ebb-dominated  system  is  characterized  by 
a  shorter  duration,  fast  ebb  current  relative  to  a  longer  duration,  slower  flood  current, 
which  is  often  associated  with  ebb-tidal  delta  morphologic  feature  (Boon  and  Byme 
1981). 

Most  of  the  original  work  on  tidal  asymmetry  has  focused  on  estuarine  systems  on 
the  U.S.  east  coast,  where  the  degree  of  the  distortion  is  primarily  described  between  the 
interaction  of  the  M2,  semidiurnal  lunar  tidal  constituent,  and  the  M4,  the  first  overtide  of 
M2  (Friedrichs  and  Aubrey  1988).  Nidzieko  (2010)  expanding  upon  the  work  of  Hoitink 
et  al.  (2003)  and  Woodworth  (2005)  demonstrated  that  the  tidal  asymmetry  for  mixed- 
mainly,  semidiurnal,  tidal  regimes  (i.e.,  short  estuaries  on  the  U.S.  west  coast)  are 
produced  primarily  through  the  combination  of  the  astronomical  lunar  (Kl),  lunisolar 
(01)  diurnal  tides  and  M2.  The  phase  relationship  between  Kl,  01,  and  M2  induces  an 
ebb-dominant  tidal  asymmetry  because  the  lower-low  water  (LLW)  follows  the  higher- 
high  water  (HHW)  (Friedrichs  1995)  and  the  phase  of  sea  surface  elevation  and  velocity 
are  in  quadrature  (Nidzieko  2010).  Though  there  is  an  ebb  dominance  at  the  estuary 
mouth,  Nidzieko  (2010)  showed  that  the  tidal  asymmetry  in  short  estuaries  is  modified  as 
a  function  of  distance  upstream  due  to  the  increased  amplitudes  of  the  overtides  and 
compound  tides,  consistent  with  the  work  by  Friedrichs  and  Aubrey  (1988).  Nidzieko 
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(2010)  provides  examples  of  estuaries,  where  the  upstream  tidal  asymmetry  (distortion) 
becomes  flood  dominated,  and  others  where  the  ebb  dominance  is  enhanced. 

An  interesting  question  is  raised  about  the  behavior  of  the  tidal  distortion  owing 
to  the  standing-wave  nature  of  a  short  estuary.  A  standing-wave  estuary  implies  that  that 
shoreward  propagating  tidal  wave  is  nearly  perfectly  reflected  at  the  landward  boundary. 
If  the  shoreward  propagating  wave  becomes  distorted  as  it  travels  landward,  does  this 
distortion  decrease  (or  increase)  as  it  reflected  back  toward  the  ocean?  To  examine  the 
shoreward  and  seaward  standing-wave  structure  of  the  tidal  wave  requires  observations 
of  both  the  potential  (pressure)  and  kinetic  (velocity)  energy.  However,  most  experiments 
tend  to  use  pressure  measurements  solely,  because  pressure  sensors  are  robust  and 
inexpensive. 

To  address  tidal  evolution  and  distortion,  a  field  experiment  was  perfonned  in  fall 
of  2011  in  Elkhorn  Slough  (ES),  CA,  which  is  one  of  the  short  estuaries  discussed  by 
Nidzieko  (2010).  Owing  to  the  standing- wave  structure,  an  array  of  four  co-located 
pressure  and  velocity  sensors  were  deployed  throughout  the  slough,  such  that  shoreward 
and  seaward  tidal  signal  can  be  separated.  The  separated  signals  allow  for  an  accurate 
description  of  the  shoreward  and  seaward  tidal  propagation  and  corresponding  distortion. 
To  date,  this  is  the  first  time  that  this  approach  is  applied  to  tidal  wave  motions  using 
field  observations.  The  amount  of  tidal  wave  reflection  for  the  astronomical  constituents, 
overtides,  and  compound  tides  are  described.  The  characteristics  of  the  shoreward  and 
seaward  tidal  wave  distortion  and  the  mechanisms  responsible  for  the  distortion  are 
discussed. 
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III.  EXPERIMENT  SETUP  AND  DATA  COLLECTION 


ES  is  located  at  the  center  of  Monterey  Bay  (MB)  and  is  adjacent  to  the  landward 
edge  of  the  Monterey  Submarine  Canyon.  The  Monterey  Bay  Sea  Floor  Mapping  Lab 
(MBSFML)  at  California  State  University  Monterey  Bay  (CSUMB)  performed  high- 
resolution  bathymetric  surveys  of  ES  between  September-December,  2011  (Figure  1).  ES 
consists  of  approximately  a  12  km  main  channel  that  is  200  m  wide  near  the  mouth  and 
narrows  upstream  to  90  m  at  Parson  Slough  (Broenkow  and  Breaker  2005).  The  length  of 
ES  is  only  3.5%  of  the  semi-diurnal  tidal  wavelength.  The  channel  from  the  edge  of  the 
man-made  jetties  to  the  HW1  bridge  is  approximately  4.5  m  deep.  East  of  the  HW1 
bridge,  there  is  an  abrupt  reduction  in  the  water  depth  to  2.5  m  (Figure  1).  There  is  no 
channel  dredging  east  of  the  bridge.  There  are  a  number  of  depth  undulations  along  the 
channel,  but  the  average  slope  is  1:4000.  The  deepest  part  of  the  channel  is  10  m  and  is 
associated  with  erosion  by  secondary  circulation  around  Seal  Bend  (Figure  1).  ES  is 
surrounded  by  extensive  Salicornia  marsh  and  mud  flats  with  a  number  of  small 
tributaries  (Broenkow  and  Breaker  2005).  At  low  tide,  the  slough  is  a  narrow  channel 
with  a  large  area  of  exposed  dry  marsh  and  mud  flats,  while  at  high  tide  the  marsh  and 
mud  flats  are  inundated. 
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Figure  1.  a)  Google  Earth  image  of  Elkhorn  Slough  showing  location  of  the 
4  ADCP  locations  with  white  dots.  The  white  dashed  line  highlights 
the  low-lying  marsh  and  mud  flats  at  the  shoreward  end  of  the  slough, 
b)  2011  bathymetry  map  of  ES  provided  by  CSUMB  SFML.  ADCP  locations 
are  shown  with  black  dots.  Depth  color-scale  provided  on  the  right. 


The  slough  is  well-mixed,  lacking  stratification  with  water  properties  closely 
resembling  MB  (Broenkow  and  Breaker  2005).  The  dominant  current  direction  in  ES  is 
streamwise  flow,  and  secondary  lateral  circulations  generated  by  curvature  in  the  bends 
of  the  slough  are  expected  to  cause  flows  to  depart  from  the  streamwise  orientation 
(Wong  1989;  Broenkow  and  Breaker  2005;  Breaker  et  al.  2008;  Nidzieko  et  el.  2008). 
Research  has  found  that  the  tidal  elevation  and  velocity  in  ES  of  the  diurnal  and 
semidiurnal  tidal  wave  are  approximately  in  quadrature  as  the  result  of  standing-wave 
characteristics  (Wong  1989;  Broenkow  and  Breaker  2005).  The  inefficient  water 
exchange  at  high  tide  due  to  the  extensive  Salicornia  marsh  and  mud  flats  and  numerous 
small  tributaries  throughout  ES  are  believed  to  enhance  longer  flood  currents  and  shorter 
more  intense  ebb  currents  (Wong  1989;  Broenkow  and  Breaker  2005).  The  wetting  and 
drying  of  the  marsh  and  mud  flats  are  hypothesized  as  the  mechanism  for  enhancing  tidal 
asymmetry  (Wong  1989;  Broenkow  and  Breaker  2005;  Nidzieko  2010). 

Four  2-MHz  Acoustic  Doppler  Current  Profilers  (ADCPs)  with  co-located 
pressure  sensors  were  deployed  longitudinally  throughout  ES  channel  for  66  days  in 
Sept-Oct  2011  (Figure  1).  The  upward-oriented  ADCPs  were  mounted  on  lm-high  Ocean 
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depth  (m) 


Science  Group  sea  spider  tripods.  The  four  sensor  locations  were  strategically  placed 
throughout  the  slough.  Station  1  and  2  (ST1,  ST2)  were  located  on  opposing  ends  of  Seal 
Bend  2.5  and  3.5  km  from  the  mouth,  respectively.  Station  3  (ST3)  was  located  upstream 
of  the  connection  to  Parson  Slough  5.5  km  from  the  mouth.  Note  that  Parson  Slough  is  an 
approximately  50  m  wide  opening  that  connects  to  a  complex  set  of  tidal  tributaries  and 
marsh  (Broenkow  and  Breaker,  2005).  Station  4  (ST4)  was  located  farther  upstream 
approximately  500  m  from  Kirby  Park  and  7  km  from  the  mouth  of  ES.  The  reflection 
point  in  ES  is  located  approximately  12  km  from  the  mouth. 

Velocities  were  obtained  at  50  cm  bin  intervals  sampled  at  1Hz  and  averaged  over 
90s  before  being  stored.  Owing  to  the  minimal  vertical  variation,  horizontal  velocities 
were  depth-averaged  and  rotated  to  the  streamwise  orientation.  The  relative  across- 
channel  variability  described  as  the  ratio  of  the  across-channel  velocity  standard 
deviation  to  the  streamwise  standard  deviation  is  8%,  13%,  7%,  and  7  %  for  ST  1-4.  The 
flow  at  ST1,  3  and  4  is  considered  streamwise  dominant,  whereas  ST2  has  more 
variability  its  proximity  to  the  secondary  circulation  associated  with  the  large  curvature 
by  Seals  Bend  (Breaker  et  al.  2008;  Nidzieko  et  al.  2008). 

The  ocean  tide  sea  surface  elevation  for  the  corresponding  ADCP  deployment 
times  was  obtained  from  the  NOAA  tidal  station  (#9413450)  located  23.5km  south  of  the 
entrance  to  Moss  Landing  in  7.3  m  of  water.  The  station  uses  a  pressure  sensor  sampling 
at  6  minutes. 
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IV.  METHODOLOGY 


To  describe  the  evolution  and  distortion  of  the  shoreward  (+)  and  the  reflected 
seaward  (-)  propagating  tidal  wave,  the  sea  surface  elevation  (r|±)  time  series  are 
constructed  from  collocated  pressure  head  (p)  and  depth-averaged  streamwise  velocity 
(u)  assuming  shallow  water, 


rf  =  \  12 


) 


(1) 


where  h  is  average  water  depth  and  g  is  the  gravitational  acceleration  (Guza  et  al.  1984). 
This  approach  successfully  described  the  cross-shore  energy  fluxes  of  infragravity  waves 
in  the  nearshore  surf  zone  environment  (Guza  et  al.  1984;  Elgar  and  Guza  1985;  List 
1992;  Sheremet  et  al.  2002).  Tidal  constituent  amplitudes  (A)  and  phases  (9)  for 
astronomical,  overtides,  and  compound  tides  are  computed  for  the  p,  and  r\±  time  series 
using  the  TTIDE  harmonic  analysis  program  (Pawlowicz  et  al.  2002).  For  a  66-day  time 
series  record,  T_TIDE  analysis  resolves  35  tidal  constituents  describing  98%  of  the 
signal,  where  2%  of  the  signal  is  considered  non-tidal  residual  motion.  Only  p  and  u  are 
used  to  separate  the  signal,  the  across-channel  velocity  (v)  is  not  considered.  Large 
variability  in  the  across-channel  velocity  (v)  can  bias  r|±  estimates.  Owing  to  the  large  v 
at  ST2  and  potential  to  bias  r|±  estimates,  results  for  ST2  will  be  omitted  unless  the 
variable  of  interest  is  p. 

In  an  effort  to  parameterize  the  tidal  evolution  and  distortion,  the  tidal  amplitudes 
are  summed  over  two  different  frequency  (f)  bands,  where  the  frequency  cut-off  between 
the  frequency  bands  is  0.1  cycles  per  hour  (cph).  T  TIDE  amplitudes  are  summed  for 
f<0.1  cph,  which  represent  the  primary  astronomical  forcing  (i.e.,  M2,  Kl,  01,  among 
others),  referred  to  as  low-frequency  motions  (subscript  lfm).  Amplitudes  summed  for 
f>0.1  cph  represent  tidal  harmonics,  such  overtides  (e.g.,  M4,  among  others),  that  are 
primarily  used  to  characterize  the  tidal  distortion,  and  are  referred  to  as  high-frequency 
motions  (subscript  hfm).  Note  that  the  total  (subscript  tot)  respresents  the  combined  sum 

of  both  the  low-  and  high-frequency  motions.  Summed  tidal  amplitudes  for  the  low- 
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frequency,  high-frequency,  and  total  motions  are  estimated  for  p  and  q±  (discussed 
below)  and  will  be  dennoted  with  subscripts. 

Tidal  amplitude  reflection  (R)  is  the  ratio  of  the  square  root  of  the  seaward  energy 
flux  divided  by  the  shoreward  energy  flux, 


Kpg(A-)-CJb  _ 

p-y 

yApg(A+)2  cgb  ]j 

py 

(2) 


where  p  is  the  density  of  sea  water,  Cg  is  the  shallow  water  wave  group  speed,  and  b  is 

the  channel  width  (Sheremet  et  al.  2002).  When  R=l,  the  tidal  wave  is  perfectly  reflected 
in  the  estuary  and  there  are  no  energy  losses  in  the  system.  An  R<  1  implies  that  there  is 
tidal  wave  energy  losses  and  R>1  implies  there  is  an  additional  source  of  energy.  Low- 
frequency,  high-frequency,  total  R  are  computed  for  ST1,  3,  and  4  (Table  1). 


ST1 

ST3 

ST4 

Rtot 

95 

102 

100 

Rffin 

90 

94 

92 

Rhfm 

125 

224 

231 

Table  1.  Reflectivity  coefficients  for  the  total,  lfm  and  hfm  for  ST1, 

ST3  and  ST4. 
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V.  RESULTS 


The  behavior  of  the  tidal  wave  is  described  by  tidal  constituent  amplitude  spectra 
for  measured  pressure  signals  at  MB  and  ST4  in  ES  (Figure  2).  Low-frequency  spectral 
amplitudes  associated  with  the  astronomical  tidal  forcing  represent  most  of  the  tidal 
variance  and  are  similar  at  both  stations.  The  high-frequency  spectral  amplitudes  are  zero 
at  MB,  but  increase  farther  upstream  in  ES.  The  high-frequency  spectral  amplitudes  are 
integer  multiple  harmonics  (i.e.,  M4,  M6)  of  the  astronomical  forcing  (i.e.,  Kl,  01,  M2) 
(Figure  2).  The  amplitude  spectra  are  consistent  with  previous  spectra  measured  ES  by 
Broenkow  and  Breaker  (2005)  and  Nidzieko  (2010). 
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Figure  2.  Tidal  constituent  amplitudes  estimated  T  TIDE  at  ST4  in  ES  from 
pressure  (black  lines)  and  the  long-term  NOAA  tidal  constituents  at  the  MB 
station.  Vertical  black  dashed  line  at  f=0.1  cph  represents  the  frequency  cut-off 
that  separates  the  low-frequency  motions  to  the  left  and  the  high-frequency 

motions  to  the  right. 


There  is  a  16%  increase  in  Atot  from  MB  to  ST4  in  ES  (Figure  3a)  estimated  by 
linear  regression.  This  trend  can  be  explained  by  the  analytical  solution  for  a  co¬ 
oscillating  tide.  Lamb  (1932)  provides  solutions  for  various  simple  geometric  inlet 
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configurations.  Of  the  configurations,  ES  is  best  described  as  constant  width  (b),  linear 
varying  water  depth  (h)  resulting  in  a  solution  defined  as, 


77  =  c 


J 0  (2-Jkx  j 
4(2 -Jkl) 


COS  (<Jt  +  s) 


(4) 


where  JQ  represents  the  zero-order,  first-kind  Bessel  function, 


,  /  is  the  slough 


length,  h0  is  the  slough  depth  at  the  end  of  the  channel,  a  is  the  radian  tidal  frequency,  x 
is  distance  from  the  landward  boundary,  C  is  the  tidal  amplitude,  s  is  the  phase 
difference.  Equation  3  applied  to  ES  geometries  results  in  an  11%  increase  in  amplitude 
from  the  mouth  to  the  head  of  the  slough,  consistent  with  the  measurements  (Figure  3a). 
While  there  is  a  difference  between  the  model  and  measurements,  the  model  describes  the 
general  co-oscillating  tidal  wave  characteristic  well  considering  the  over  simplification  of 
the  slough  (constant  width). 


0  1000  2000  3000  4000  5000  6000  7000  8000 

x  (m) 


Figure  3.  a)  Squares  represent  Atot  between  MB  and  ST4  as  a  function  of 
along-channel  distance.  The  dashed  line  is  the  best-fit  line  for  Atot  estimated  by 
linear  regression.  Solid  line  represents  the  Lamb  (1932)  analytical  solution  for 
co-oscillating  tide  for  ES.  b)  Ap,hlm  (squares)  between  MB  and  ST4  as  a  function 
of  along-channel  distance.  The  dashed  line  is  the  best-fit  line  for  Ap,hfm  estimated 
by  linear  regression,  c)  Rtot  (triangles),  Rum  (circles),  and  Rhfm  (squares)  as  a 
function  of  along-channel  distance. 
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The  Ap  hfm  (estimated  with  pressure)  steadily  increases  from  0  to  approximately 
14  cm  between  MB  and  ST4  (Figure  3b).  The  increase  in  the  tidal  wave  non-linearity 
with  along-channel  distance  associated  with  pressure  signal  is  consistent  with  previous 
research  (Speer  and  Aubrey  1985;  among  others).  However,  the  measured  pressure  signal 
represents  the  potential  energy  portion  of  a  tidal  wave.  A  purely  progressive  wave  can  be 
described  solely  by  pressure  or  velocity  because  they  are  in  phase.  For  a  standing-wave, 
both  potential  and  kinetic  energy  must  be  considered  to  effectively  characterize  the  wave 
(i.e.,  Equation  1). 

The  reflection  coefficients,  R  (Equation  2),  for  each  station  for  the  total,  low- 
frequency,  and  high-frequency  related  shoreward  and  seaward  tidal  signal  provide 
insights  to  the  bulk  tidal  wave  characteristics  and  its  evolution.  The  Rtot  ranges  from 
95%-102%  (Figure  3c,  Table  1).  The  small  differences  in  Rtot  are  related  to  the  fact  that 
the  propagating  shoreward  and  seaward  tidal  can  behave  different  at  the  different 
stations.  Owing  to  this,  ST1  represents  the  overall  tidal  behavior  in  the  slough,  because  it 
is  located  the  farthest  from  the  upstream  boundary  and  describes  the  reflection  over  a 
longer  distance,  whereas  ST3-4  represents  the  tidal  behavior  farther  upstream.  The  Rtot 
results  imply  that  farther  upstream  the  tide  is  more  reflected  than  downstream.  For  ST1, 
5%  of  the  total  amplitude  is  dissipated  presumably  by  bottom  friction  (Speer  and  Aubrey 
1985)  or  secondary  circulation  associated  with  Seal  Bend  (Broenkow  and  Breaker  2005; 
Breaker  et  al.  2008;  Nidzieko  et  al.  2008).  In  addition,  the  marsh  and  mud  flats  of  Parson 
Slough  located  between  ST1  and  ST3  may  contribute  to  the  reflection  at  ST1  and  its 
difference  measured  at  ST3-4.  ST3  and  ST4  are  closer  to  the  reflection  point  and  the 
complexities  of  the  slough  are  much  less,  as  there  are  minimal  tidal  tributaries  and  mostly 
marsh  and  mud  flats.  The  high  Rtot  at  each  station  suggest  ES  is  highly  reflective, 
consistent  with  Broenkow  and  Breaker  (2005),  Breaker  et  al.  (2008)  and  Nidzieko 
(2010). 

Rifmat  ST1,  3  and  4  are  90%,  94%  and  92%,  while  Rhfm  at  ST1,  3  and  4  are  125%, 
224%  and  231%,  respectively  (Figure  3c,  Table  1).  Rifm  indicate  approximately  10% 
energy  loss  compared  to  the  gain  in  energy  associated  with  Rhfm.  Owing  to  the  fact  that 
Rtot  represents  contributions  from  Rifm  and  Rhfm,  and  that  Rtot  is  -100%,  then  the  A,rjllm  of 


13 


the  seaward  wave  must  be  more  non-linear  than  the  An_  hfm  of  the  shoreward  wave.  The 
larger  A11+,hfm  of  the  seaward  wave  results  in  a  Rhfm  greater  than  1 00%  and  accounts  for 
the  10%  energy  loss  associated  with  Rifm.  The  larger  Rhfm  value  suggests  energy  is 
actually  transferred  from  the  low-frequency  tidal  motions  (i.e.,  astronomical  constituents) 
to  the  high-frequency  tidal  motions  (i.e.,  harmonics)  (Parker  1991).  These  are  new 
findings  previously  unavailable  from  only  pressure  measurements. 
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VI.  DISCUSSION 


The  previous  section  highlighted  the  differences  between  the  shoreward  and 
seaward  waves  using  bulk  statistics,  but  the  bulk  statistics  do  not  describe  and  explain 
why  the  waves  are  actually  different.  Many  researchers  believe  that  tidal  distortion  in 
estuaries  with  a  main  channel  and  vast  low  lying  marsh  and  mud  flats  is  due  to  inefficient 
water  exchange  around  high  tide  as  the  marsh  and  mud  flats  become  inundated 
(Nummedal  and  Humphries  1978;  Boon  and  Byrne  1981;  Speer  and  Aubrey  1985; 
Friedrichs  and  Aubrey  1988  Wong  1989;  Broenkow  and  Breaker  2005;  Breaker  et  al. 
2008;  Nidzieko  2010).  Unfortunately,  field  measurements  in  this  experiment  were  not 
obtained  at  the  landward  boundary  to  describe  the  effects  of  the  marsh,  but  the  time 
signals  of  the  shoreward  and  seaward  tidal  wave  are  available  to  provide  insight  on  the 
behavior  that  occurs  at  the  landward  boundary,  which  is  examined  next. 

The  characteristics  of  a  measured  24-hr  time  signal  for  yearday  269  of  the 
shoreward  (Figure  4,  solid  lines)  and  seaward  (Figure  4,  dashed  lines)  tidal  wave  for  ST1, 
ST3  and  ST4  are  evaluated.  The  shoreward  tidal  wave  first  propagates  by  ST1  (Figure  4, 
blue  line)  and  then  to  ST3  (Figure  4,  green  line)  and  onto  ST4  (Figure  4,  red  line),  which 
is  shown  by  the  temporal  offset  between  the  stations.  The  temporal  lag  between  each 
station  represents  the  propagation  time  that  is  dependent  upon  the  distance  between  the 
stations  (e.g.,  the  largest  distance  is  between  ST1  and  ST4  resulting  in  the  largest 
temporal  difference).  The  seaward  time  signals  (Figure  4,  solid  lines)  are  similar  to 
shoreward  time  signals  (Figure  4,  dashed  lines),  except  in  reverse.  The  seaward  time 
signal  represents  the  tidal  wave  that  has  propagated  to  the  landward  boundary  and 
reflected  downstream.  Therefore,  ST4  represents  the  first  station  that  the  tidal  seaward 
wave  encounters,  and  then  the  tidal  wave  propagates  past  ST3  and  then  ST1.  The  time 
difference  between  the  shoreward  signal  and  the  seaward  signal  represents  the  time  the 
tidal  takes  to  propagate  to  the  landward  boundary  and  back,  which  is  referred  to  as  the 
reflection  time.  The  reflection  time  increases  with  increasing  distance  downstream  from 
the  landward  boundary.  The  reflection  time  is  largest  for  ST1. 
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Figure  4.  The  shoreward  (solid  lines)  and  seaward  (dashed  lines)  sea  surface 
elevation  measured  for  yearday  269  (shown  in  relative  hours)  for  ST1  (blue), 

ST3  (green),  and  ST4  (red). 

The  time  difference  between  the  shoreward  and  seaward  time  signal  steadily 
increases  (Figure  4)  from  low  tide  to  high  tide  implying  that  the  reflection  time  is 
changing  and  appears  to  be  a  function  of  elevation  (total  depth).  The  reflection  time 
decreases  from  high  tide  to  low  tide  and  the  response  differs  slightly  from  station  to 
station  (Figure  4).  The  reflection  time  for  high  tide  maxima  and  low  tide  minima  are 
estimated  for  the  66-day  record  for  ST1,  ST3,  and  ST4,  resulting  in  four  estimates  per 
day.  A  linear  trend  is  found  between  the  reflection  time  and  the  sea  surface  elevation, 
which  increases  with  increasing  elevation  (Figure  5).  The  linear  slopes  are  0.83,  0.99,  and 
0.87  for  ST1,  ST3,  and  ST4  suggesting  about  a  1  hr  time  lag  per  meter  of  elevation  gain. 
Though  a  linear  trend  is  observed,  there  is  a  variability  indicated  by  r  -values  of  0.59, 
0.69,  and  0.60  for  ST1,  ST3,  and  ST4.  The  large  scatter  at  ST1  (Figure  5a)  may  be  related 
to  influence  of  Parson  Slough  or  Seal  Bend.  ST3  and  ST4  interestingly  have  reduced 
variability  at  low  tide  (Figure  5b, c)  suggesting  that  the  reflection  time  below  MSL  is 
controlled  by  the  channel  and  above  MSL  controlled  by  the  more  complex  nature  of  the 
mud  flats  and  marshes. 
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Figure  5.  The  reflective  time  lag  estimated  from  daily  high  tide  and  low  tide 
maxima  and  minima  for  a)  ST1,  b)  ST3,  and  c)  ST4.  Dashed  line  represents 
the  best-fit  line  estimated  with  linear  regression.  The  linear  slopes  are  0.83,  0.99, 
and  0.87  and  the  revalues  are  0.59,  0.69,  and  0.60  for  ST1,  ST3,  and  ST4. 

These  results  are  initially  counter-intuitive,  because  it  is  expected  as  the  depth  of 
the  main  channel  increases  (decreases)  with  a  rising  (falling)  tide,  the  phase  speed, 

yjg(h  +  rj)  ,  of  the  tidal  wave  would  increase  (decrease)  [Friedrichs  and  Aubrey  1988; 
Wong  1989;  Broenkow  and  Breaker  2005]  and  decrease  the  reflection  time  which  differs 
with  the  observations.  The  total  water  depth  (h+p)  assumes  that  the  channel  has  vertical 
walls.  While  this  is  relatively  true  for  some  portions  of  the  slough  located  downstream,  at 
the  upper  end  of  the  slough  the  channel  is  surrounded  by  shallow  marsh  and  mud  flats 
(Figure  1),  which  becomes  inundated  at  higher  tides  (Broenkow  and  Breaker  2005).  Since 
the  channel  cannot  be  described  by  simple  rectangle,  the  average  (effective)  water  depth 
is  estimated  by 


where  Ac  is  the  cross-sectional  area  and  b  is  the  channel  width.  Therefore,  increasing  tidal 

elevation  results  in  a  decrease  in  average  water  depth  as  water  spreads  out  over  the 
marsh.  The  result  is  a  decrease  in  phase  speed  and  increased  reflection  time  as  observed 
(Figures  4  and  5).  This  is  consistent  with  the  results  found  for  a  progressive  wave  system 

by  Friedrichs  and  Aubrey  (1988),  where  the  volume  of  intertidal  storage  divided  by  the 
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channel  volume  is  large  in  relation  to  the  tidal  amplitude  divided  by  water  depth  that 
caused  the  wave  at  high  tide  to  propagate  slower  than  at  low  tide. 

The  difference  in  reflection  time  as  function  of  elevation  modifies  the  seaward 
tidal  signal  relative  to  the  shoreward  time  signal,  which  explains  why  the  seaward  tidal 
wave  becomes  more  distorted.  If  the  reflection  time  is  constant,  then  the  shoreward  and 
seaward  waves  would  be  the  same.  However,  this  cannot  happen  even  in  a  rectangular 
channel  because  the  phase  speed  would  increase  with  increasing  depth  (Friedrichs  and 
Aubrey  1988).  This  suggests  that  the  seaward  tidal  wave  will  always  be  more  distorted 
than  the  shoreward  tidal  wave  in  a  reflective  estuary. 
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VII.  SUMMARY  AND  CONCLUSIONS 


Co-located  measurements  of  pressure  and  depth-averaged  velocity  measured 
longitudinally  throughout  a  short  slough,  allowed  the  tidal  signal  to  be  separated  into 
shoreward  and  seaward  tidal  signals,  providing  new  insights  into  tidal  wave 
characteristics.  The  overall  tidal  wave  reflection  in  Elkhorn  Slough,  CA,  is  high  (-100%). 
The  low-frequency  seaward  tidal  signal,  associated  with  astronomical  tidal  constituents, 
tends  to  lose  10%  of  its  amplitude  relatively  to  the  shoreward  tidal  signal,  whereas,  the 
high-frequency  seaward  tidal  signal,  associated  with  tidal  harmonics,  has  a  large  (>125%) 
increase  in  amplitude  relative  to  the  shoreward  tidal  signal.  The  increase  in  amplitude  of 
the  high  frequency  seaward  tidal  signal  suggests  energy  is  transferred  from  the 
astronomical  tidal  constituents  to  the  harmonics  and  overtides  consistent  with  the 
findings  of  Parker  (1991).  The  high-frequency  amplitudes  are  responsible  for  the  tidal 
wave  distortion  resulting  in  the  seaward  tidal  wave  being  more  distorted  than  the 
shoreward  tidal  wave.  The  mechanism  responsible  for  the  seaward  tidal  wave  distortion 
can  be  explained  by  the  shallow-water  wave  phase  speed,  which  is  depth  dependent.  As 
the  tidal  elevation  increases,  the  small  channel  near  the  landward  boundary  fills  up 
quickly  and  water  flows  onto  the  adjacent  shallow  mud  flats  and  marshes  reducing  the 
effective  depth  of  the  channel.  The  decreases  in  the  effective  depth  reduces  the  shallow- 
water  wave  phase  speed  inducing  a  reflective  time  lag,  defined  as  the  time  the  tidal  wave 
takes  to  propagate  upstream  to  the  landward  boundary  and  back.  A  near  linear  trend  was 
found  for  the  reflective  time  lag  that  increases  with  increasing  elevation  for  all  stations. 
This  trend  effectively  modifies  the  characteristic  of  the  seaward  tidal  signal  relative  to  the 
shoreward  tidal  signal,  which  causes  the  tidal  wave  distortion.  In  summary,  the  reflective 
behavior  of  the  slough  is  responsible  for  the  increase  in  tidal  distortion  for  the  seaward 
tidal  wave. 
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